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ncutralirio production, measurable in this channel. The + £ signature is investigated, where 
i is an electron or a muon, with the additional requirement of large missing transverse energy. In 
this analysis, the leptori transverse momenta with respect to the beam direction (pr) are as low 
as 5 GeV/c, a selection that improves the sensitivity to particles which are light as well as to ones 
which result in leptonically decaying tau leptons. At the same time, this Iow-jjt selection presents 
additional challenges due to the non-negligible heavy-quark background at low lepton momenta. 
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and trilepton control regions are investigated, and good agreement between experimental results and 
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I. INTRODUCTION 

The standard model (SM) of particle physics is enor- 
mously successful in describing the known particles and 
their interactions. However, strong motivation from ex- 
perimental data as well as important theoretical consid- 
erations point to new physics beyond the SM. Astrophys- 
ical observations that have resulted in the "concordance" 
model of cosmology [l| require a source of dark matter 
that does not exist in the SM. Theoretically, the SM has 
well-known limitations in explaining the origin of mass 
and solving the hierarchy problem. Moreover, it does 
not satisfy our desire for the unification of the strong and 
electroweak interactions and the integration of gravity in 
a unique theory 0. 

A powerful strategy for discovering new physics is to 
search in event topologies where the SM predicts ex- 
tremely low production rates. One of these topologies 
is three leptons (trileptons) in hadronic collisions. The 
lepton candidates we observe at the Tevatron collider re- 
sult mainly from QCD processes or the decays of massive 
gauge bosons {W or Z), or photon conversions. The lep- 
tons rarely appear in multiplicity greater than two. The 
trilepton signature is favored by a large class of models 
of supersymmety (SUSY) [1, 0], in which the hghtest su- 
persymmetric particles are the gauginos, the supersym- 
metric partners of the gauge bosons. The corresponding 
observable SUSY particles are two charginos {xt2) ^^'^ 
four neutralinos (x? 2 3 4)1 which result from the mixing 
of the gauginos and the supersymmetric partners of the 
Higgs bosons, the higgsinos. The associated production 
of charginos and neutralinos may have a detectable cross 
section Q at the Tevatron and may give rise to trilepton 
events as shown in Figs. [T] and [H The most common de- 
cays are through off-shell vector bosons or scalar leptons 
(sleptons), with branching fractions that depend on the 




q Xi 

AAAAAAAAy" 

q' xi 

*\AAAAAAA/ 

FIG. 1: Chargino-neutralino production through an s-channel 
W boson (a) and a t-channel squark propagator (b). The t- 
channel is suppressed in scenarios with very massive squarks. 




FIG. 2: Chargino and neutralino decays through gauge bosons 
(a,b) and through sleptons (c,d). The branching fractions 
depend on the masses of the sleptons, which always decay 
to charged leptons, unlike the gauge bosons. The leptonic 
signature consists of three leptons in both cases. 
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chargino, neutralino, and slepton masses. 

Under the assumption of i?-parity [6| conservation, 
SUSY particles cannot yield only SM particles in their 
decay; the lightest SUSY particle (LSP) will be stable 
and escape detection. Therefore, SUSY events would be 
characterized by large transverse momentum imbalance 
("missing transverse energy", or I^t)- In many SUSY 
scenarios the lightest neutralino is either the LSP or it de- 
cays to the LSP resulting in in both cases. Additional 
l^T results from the undetected final-state neutrinos, as 
shown in Fig. [2l The trilepton-|-^T topology investigated 
here is the "golden" signature for the discovery of SUSY 
at the Tevatron 0, S, S [M [H ■ 

The LSP is a candidate for the cold dark matter of the 
universe [H, • In addition SUSY offers a solution to 
the hierarchy problem [3, [l^ 13 ^^'^ ^be possibility for 
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unification of interactions at high energies [l7l |. 

Searches for chargino and neutraUno production have 
been previously performed by the LEP [H, [l^ and Teva- 
tron experiments [20L [2l| . In this paper we present 
a trilepton analysis that utilizes increased luminosity 
and improved kinematic acceptance. We search for new 
physics in the final state with two muons and an ad- 
ditional electron or muon using data collected with the 
CDF II detector from March 2002 to February 2006 from 
proton- antiproton collisions at y/s = 1.96 TeV. The inte- 
grated luminosity of our sample is 976 pb~^. To increase 
the sensitivity to new light particles and tau leptons that 
decay leptonically, we use a very low pT threshold (5 
GeV/c) for the identified leptons, where pt is the trans- 
verse momentum with respect to the beam direction. 

We define several dimuon and trilepton SM-dominated 
control regions, in which we verify our understanding of 
the backgrounds. In order to avoid bias, we complete 
the validation of the background - both in event yields 
and kinematic shapes - in the control regions before in- 
vestigating the events in the signal region. Finally, our 
result is combined with other trilepton searches at CDF 
[2^ to set a stronger limit on chargino-neutralino produc- 
tion. Although this search is inspired by SUSY-predicted 
chargino-neutralino production, the analysis is generic 
enough to be sensitive to any new physics that would 
enhance the production of prompt trileptons and I^t- 

This paper is organized as follows. In Section[ll]we de- 
scribe the CDF II detector. In Section ITTTl we define the 
experimental dataset and present an event selection that 
reduces the SM background expectation while accepting 
events from possible new-physics signals. Section |IV] de- 
scribes how the SM background rates are estimated, and 
Section |V] discusses two SUSY-model scenarios we con- 
sider. In Section IVTl we present the determination of the 
systematic uncertainties on the signal and background 
event-yield predictions. In Section IVIII we present the 
event yields and kinematic distributions in our control 
regions that increase the confidence in our understand- 
ing of the SM background. Finally, in Section I Villi we 
present the results in the signal region. 



II. THE CDF II DETECTOR 

The CDF II detector [2^ is a multi-purpose cylindri- 
cal detector with projective-tower calorimeter geometry 
and excellent lepton identification capability. It oper- 
ates at the Tevatron collider where protons and antipro- 
tons collide with a center-of-mass energy of 1.96 TeV. In 
our coordinate system the positive z-axis is defined by 
the proton beam direction and the positive y-axis by the 
vertical upward direction. The detector is approximately 
symmetric in the rj and (j) coordinates, where the pseudo- 
rapidity 77 is defined as 77 = — ln(tan(6'/2)), is the polar 
angle with respect to z, and (p is the azimuthal angle. We 
briefiy present here the CDF components that are most 
critical to this analysis. 



In the center of the apparatus, near the beam collision 
point, a silicon detector of inner radius of 1.35 cm and 
outer radius of 25.6 cm provides detailed tracking in the 
I77I < 2 region, necessary for the accurate determination 
of the proton-antiproton interaction points (primary ver- 
tices) and impact parameters of particle trajectories with 
respect to these points. 

A cylindrical 96-layer open-cell argon-ethane (50%- 
50%) drift chamber (COT) of inner radius of 44 cm and 
outer radius of 132 cm provides tracking for charged 
particles with ~ 100% detection efficiency in the cen- 
tral (jryl < 1.1) region. The central tracking system 
is located in a magnetic field of 1.4 T provided by a 
superconducting solenoidal magnet. The relative reso- 
lution in tracking momentum provided by the COT is 
Spt/pt = 0.0017]3T(GeV/c)-i. 

Surrounding the central tracker, and outside the 
solenoid, a central electromagnetic calorimeter (CEM) 
and a central hadronic calorimeter (CHA) measure the 
energy of electrons, photons, and hadrons. The CEM 
is composed of layers of lead and scintillator whereas 
the CHA is composed of layers of steel and scintilla- 
tor. The relative energy resolution is 13.5%/-\/Et © 2% 
for the CEM and 75%//E^ ® 3% for the CHA, where 
the transverse energy Et = EsmO is quoted in GeV 
units. A strip chamber (CES), placed inside the elec- 
tromagnetic calorimeter at the position of maximum de- 
velopment of the electromagnetic shower (six radiation 
lengths), is used for shower shape determination and for 
matching the calorimeter energy depositions with COT 
tracks. In the forward region, a plug electromagnetic 
calorimeter (1.1 < \t]\ < 2.4) has a relative resolution 
of 16%/a/Et © 0.7% and a plug hadronic calorimeter 
(1.3 < \t]\ < 2.4) a resolution of UO%/^/E^ ® 4%. 
The raw missing transverse energy vector is defined as 
— (X^i-E-jr), where Eix has magnitude equal to the en- 
ergy deposited in the i"^ calorimeter tower and direc- 
tion perpendicular to the beam axis and pointing to that 
calorimeter tower. 

Outside the calorimeters, the central muon system 
consists of drift chambers. The central muon cham- 
bers (CMU) detect muons in the pseudorapidity range 
\ri\ < 0.6, while the central muon extension (CMX) cham- 
bers detect muons in the 0.6 < \ri\ < 1.0 range, both with 
a detection efficiency of almost 100% for muons above 3 
GeV/c. To reduce the hadron punch-through contami- 
nation, extra chambers (CMP) are installed outside the 
CMU chambers, with extra steel absorber added between 
them. The muons that are detected by both CMU and 
CMP chambers are labeled "CMUP muons" , and their 
detector signatures cannot be easily caused by hadrons. 

The instantaneous luminosity is measured with 
Cherenkov counters located close to the beam line at 
3.7 < |?7| < 4.7. 

The CDF trigger system [23| has a three-level pipelined 
and buffered architecture; each level provides a rate re- 
duction sufficient to allow for processing at the next level 
with minimal deadtime. The first level consists of special- 
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purpose processors that accept events at rate of 25 kHz, 
with an average event size of 170 kB, counts main trig- 
gering objects and feeds the second level with an event 
rate of 350 Hz. The second level is also based on hard- 
ware and performs a partial event reconstruction before 
passing the events to the next level. Finally a software- 
based third level uses a fast version of the offline event 
reconstruction to reduce the event rate to 75 Hz, ap- 
propriate for writing to tape. The track-based triggers 
account for approximately 75% of the trigger bandwidth 
and are used in this analysis. For a muon trigger, the 
main requirement is that a COT track is geometrically 
matched to a track segment in a muon detector. 



III. THE CDF DATASET AND 
SIGNAL-REGION EVENT SELECTION 

In order to include in our analysis muons and elec- 
trons that come from tau decays, we use a low trans- 
verse momentum requirement {pt > 5 GeV/c) for these 
leptons. For this reason we analyze data collected with 
the CDF low-pT diniuon triggers {pr nominally above 
4 GeV/c for both muons). These muons are central in 
the detector (CMUP or CMX). We measure the trig- 
ger efficiency using J/ip, T and Z-boson events collected 
with single-muon triggers. In these samples, we remove 
hadronic backgrounds using the mass-spectra sidebands, 
and count the frequency that a second muon fired the 
trigger of interest. The plateau value of the trigger effi- 
ciency's pt dependence for single muons is ^ 0.95 and it 
is reached a.t px ^ 5 GeV/c. 

After the collected events are processed by the offline 
reconstruction software, additional requirements are ap- 
plied for the definition of the dimuon sample. We re- 
quire that each event has a primary vertex within 60 cm 
from the nominal center of the detector in the z direction 
and that at least two muons with transverse momenta 
above 5 GeV/c originate from that primary vertex and 
pass the CDF standard muon tracking and calorimetry 
requirements and track-chamber matching requirements 
[25|. In events with more than one reconstructed pri- 
mary vertex, we use the primary vertex that is closest to 
the tracks of the two highest-pT muons that satisfy all 
other event requirements. We specifically require that 
two good-quality COT tracks are geometrically matched 
with respective reconstructed track segments in the CMX 
or CMU-I-CMP detectors, that the energies deposited in 
the electromagnetic and hadronic calorimeters are con- 
sistent with that expected from minimum ionizing par- 
ticles, and that the two muons are isolated. We define 
the isolation / as the energy deposited in the calorime- 
ters in a cone of AR = a/ (Ai/))^ -t- (Ar;)^ = 0.4 around 
the muon without counting the energy deposited by the 
muon. We require that J < 0.1 x ptc if pt > 20 GeV/c 
or 7 < 2 GeV otherwise, where pt is the transverse mo- 
mentum of the muon. The selected two muons are also 
AR > 0.4 apart. A critical requirement is that the muons 



are prompt as measured by the impact parameter {do), 
defined as the distance of closest approach of a track to 
the primary vertex in the transverse plane. We require 
that |do| < 0-02 cm if the the muon leaves tracking signals 
in the silicon detector (silicon hits) and that jdo] < 0.2 cm 
if the muon leaves no silicon hits. We expect that most 
muons with large impact parameters come from heavy 
flavor (bottom- or charm-hadron semileptonic decays), 
fake muons (light-flavor hadrons such as pions and kaons 
that decay in flight or punch-through to the muon de- 
tectors), and cosmic rays. The heavy flavor (HF) and 
fake- muon backgrounds dominate at low dimuon masses. 
Residual cosmic-ray background, not removed by the cos- 
mic filters described in [2^, is reduced by requiring that 
the three-dimensional angular separation (Aip) of the two 
highest-pT muons is less than 178 degrees. After includ- 
ing the selection criteria discussed above, the total muon 
identification efficiency, as measured with J/ip and Z bo- 
son CDF data, is (90-96)%, rising with increasing muon 
Pt- 

For the trilepton selection, we require the presence 
of a third muon satisfying the same selection require- 
ments as the first two, or an electron satisfying the 
CDF standard electron calorimeter, tracking, and track- 
calorimeter matching identification requirements [2^ . 
The transverse energy and momentum of an electron is 
required to exceed 5 GeV. Tracks associated with elec- 
trons should match hits in the CES wires. We require 
that / < 0.1 X Et if Et > 20 GeV or / < 2 GeV oth- 
erwise, where / is now the energy-based isolation of the 
electron, and Et is its transverse energy. Electrons orig- 
inating from photons that convert into e+e~ pairs are 
identified with an algorithm [25] that seeks nearby tracks 
with a common vertex and direction. These electrons 
are removed from the observed data sample. The elec- 
tron identification efficiency is (75-83)% [25, rising with 
increasing electron transverse energy, as measured with 
Drell-Yan [2^ electrons. The third lepton is required to 
be AR > 0.4 away from the leading two muons. 

We define the signal region by the following additional 
requirements: the dimuon mass (constructed using the 
two highest-pT muons) is greater than 15 GeV/c^, for 
removal of low-mass resonances, and outside a Z mass 
window of 76 < M^^ < 106 GeV/c^. In addition, we 
require the missing transverse energy (^t) to exceed 15 
GeV, in order to select events with undetected new parti- 
cles while rejecting Drell-Yan, HF, and fake-muon back- 
grounds. Finally, we count the number of jets A'jets with 
energy above 15 GeV and we require that A'jcts < 1, in 
order to reduce the ti background. In this analysis we 
use jets defined by a fixed-cone algorithm [2^ with a 
cone size of AR — 0.4. We require that jets deposit less 
than 90% of their measured energy in the electromag- 
netic calorimeter, in order to avoid countin g el ectrons or 
photons as jets. Jet energies are corrected [27[ to repre- 
sent better the energy of the final-state hadrons. Global 
and local corrections are applied as well as inclusion of 
corrections for the effects of multiple interactions. These 
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corrections are also applied to the raw missing transverse 
energy for the calculation of I^t, which is also corrected 
for the presence of muons in our events. We check the 
consistency of the observed data compared to the SM 
predictions in the control regions that are described in 
Section IVni 



B. Data-estimated backgrounds 

We first estimate the fake-lepton background, using 
an independent CDF data sample. Subsequently, we use 
this fake-lepton background and the MC-estimated DY 
contribution in our HF-estimation method. 



IV. STANDARD-MODEL BACKGROUNDS 

To determine the significance of any incompatibility 
between prediction and observation, and also to set lim- 
its on production cross sections and masses of new parti- 
cles, we need a reliable background estimation. The ma- 
jor SM source of dimuons is the Drell-Yan (DY) process 
and, in events with low dimuon mass, HF production and 
the fake-lepton background. In the trilepton regions, the 
dominant backgrounds are DY (accompanied by a fake 
lepton), dibosons {WW, ZZ, and WZ), and HF. Because 
HF and fake leptons are difhcult to model with Monte 
Carlo (MC) simulations due to sizable higher-order QCD 
effects and the imperfect modeling of the lepton isolation 
in a high particle- multiplicity hadronic environment, we 
estimate these backgrounds using CDF data. All other 
backgrounds are estimated with MC simulation. 



A. MC-estimated backgrounds 

We use the pythia [28i] generator to model the DY, 
WW, ZZ, and tt background, and madevent [2^ for 
the WZ background [301 ■ The DY background includes 
the decays to tau leptons that subsequently decay to 
muons [Slj. We use the CTEQSl [s^I parton distribution 
fimctions (PDF) throughout. For the trilepton predic- 
tions we require the reconstructed electrons and muons 
to be kinematically matched with the generator-level lep- 
tons, in order not to double-count some of the fake-lepton 
contribution. To estimate the trilepton background from 
DY-I-7, we relax this matching requirement, demand that 
the electron is identified at the event-simulation level as a 
photon-conversion product, and normalize the surviving 
event using a scale factor 33] . This scale factor accounts 
for the difference in conversion-removal inefficiency be- 
tween the observed data and the MC simulation. In the 
remainder of the paper we add the DY-I-7 background 
to the rest of the diboson contribution {WW, ZZ, and 
WZ). We process each generated event with the CDF 
detector simulation, based on GEANT [s^I- We normal- 
ize all samples using the leading-order theoretical cross 
sections multiplied by the appropriate scaleO'K-factor" ) 
to correct for next-to- leading order effects [H, . Scale 
factors that correct for the known differences in lepton 
identification and reconstruction efficiencies between the 
observed data and the MC simulation are also applied. 



1. Fake leptons 

"Fake" leptons are reconstructed lepton candidates 
that are either not real leptons or are real leptons but are 
neither prompt nor do they originate from semileptonic 
decays of HF quarks. In the case of muons, the fakes can 
be light-flavored hadrons, such as pions and kaons or part 
of hadronic showers, that penetrate ("punch through") 
the calorimeters and reach the muon detectors or decay 
to muons in flight. In the case of electrons, fakes are 
jets that are misrcconstructcd as electrons, often due to 
neutral pions that decay to photons, which shower in 
the electromagnetic calorimeter. We can thus associate 
the fake leptons with light-flavor partons. Using multijet 
CDF datasets collected with jet-based triggers, we mea- 
sure the "fake rate", i.e., the probability for an isolated 
track to be misreconstructed as a muon or the probabil- 
ity for a jet to be misreconstructed as an electron. The 
fake-rate is measured as a function of the track's (jet's) 
transverse momentum (energy) and pseudorapidity. The 
fake rate is of the order of 10~^ for isolated tracks to be 
incorrectly reconstructed as muons and increases with the 
Pt of the muon candidate's track. The fake rate of a jet 
being reconstructed as an electron is of the order of 10"''' 
and falls with increasing Et. The fake rates increase for 
higher pseudorapidity leptons [s^. 

To determine the background coming from a real 
muon and a misidentified hadron (i.e., fake dimuon back- 
ground), we use single- muon low-pT-triggered CDF data. 
For each event, we require one good muon candidate that 
passes the requirements of our analysis. We then apply 
the fake rate on all other tracks in the event, except on 
the track of a second muon (to remove DY contamina- 
tion of the fake background). We remove events in which 
the "muon-|-track" mass is within the Z boson window 
(76 < Mf,f, < 106 GeV/c^) and also 35 < < 55 
GeV. These events are associated with decays of real Z 
bosons produced at rest, where one decay muon is not 
detected, resulting in equal to about half the mass of 
the Z boson. We investigate the heavy flavor contami- 
nation in this light-flavor-dominated background, caused 
by a real muon coming from a heavy-quark semileptonic 
decay, which is misreconstructed as an isolated track in- 
stead of a muon. This contamination is negligible (ap- 
proximately 0.2% of the background), mainly because of 
the background-rejection power of our muon isolation re- 
quirement. Because the single-muon low-px trigger was 
not always present during data taking, the dimuon fakes 
extracted using this trigger are normalized to the default 
dimuon CDF data luminosity, in order to represent the 



8 



size of fake dimuon contamination in our analysis CDF 
dataset. In CDF data with no Ipx or jet multiplicity cuts 
applied ("inclusive" dataset), ~ (9 ±5)% of the dimuons 
are fake. In the signal region, the dimuon fake contami- 
nation is ~ (16 ± 8)%. 

For the determination of the background coming 
from a real muon pair and a misidentified hadron (i.e., 
fake trilepton background), we use our dimuon low-p^- 
triggered CDF dataset, require two good muons, and 
model the fake third lepton by applying the fake rate to 
the extra tracks. We assume that the number of events 
with two fake leptons is negligible, given the low value of 
the fake rates. In order not to over-count the trilepton 
fakes, we require the three leptons in our signal MC and 
background MC samples to be kinematically matched 
with the generated ones. The fake trilepton background 
is determined to be ^ (50±25)% of the total background 
in both the inclusive dataset and the signal region. 



2. Heavy flavor 

One of the most significant challenges of this analysis 
is the consideration of muons with transverse momentum 
as low as 5 GeV/c. This low pT requirement increases 
our acceptance, but at the same time contaminates our 
sample with HF and fake-lepton events. 

We present here an innovative technique for the de- 
termination of the amount of this background using the 
observed data. We construct an HF-rich (HFR) CDF 
dataset by reversing the impact parameter requirement 
for at least one of the observed muons, so that the abso- 
lute value of the muon impact parameter is above 0.02 cm 
if there are silicon hits associated with the muon track, or 
above 0.2 cm if there are no silicon hits. We also require 
the dimuon mass to be less than 35 GeV/c^. Monte Carlo 
studies show that above that value we expect mainly DY 
and a negligible HF background. We investigated the ex- 
pected dimuon mass spectrum of DY and fake-lepton in 
the HFR sample and we determined that the effect of the 
contamination is negligible. 

We subsequently use the HFR dimuon mass shape 
combined with the absolute fake dimuon mass distribu- 
tion plus the absolute DY dimuon mass distribution from 
MC simulation in order to fit the observed data. All 
data samples other than HFR include the low impact 
parameter requirement. Because we observe negligible 
DY in the same-charge dimuon channel, we perform the 
fit for same-charge and opposite-charge dimuons sepa- 
rately. This helps us validate our HF-estimation method 
in the HF-rich same-charge dimuon environment. The 
only free parameter of the fits is the HF normalization 
- the DY contribution is fixed based on the theoretical 
cross section and the integrated luminosity of the ob- 
served data, and the fake-lepton background distribution 
is fixed based on the absolute expectation, as described in 
Section llVB II The results of the fits can be seen in Figs. 
[3]and|31 for opposite-charge and same-charge muons, re- 



spectively. From the two fits, we extract two HF normal- 
ization factors that are applied as weights to the origi- 
nal unweighted same-charge and opposite-charge dimuon 
HFR events, in order to describe the HF background in 
the observed data. The same weights are used in all 
kinematic control regions. The weight for the opposite- 
charge HF dimuons is 1.94±0.04 and for the same-charge 
HF dimuons is 1.12 ± 0.05, where the uncertainties come 
from the fits. This tells us that we expect almost twice as 
many opposite-charge HF events in the low-impact pa- 
rameter region compared to the high-impact parameter 
one. Overall, the ratio of opposite-charge HF to same- 
charge HF events in the inclusive observed data after 
the normalization is ^ 4:1, a value that is also verified 
with bb/ cc MC simulation and is a result of the conserved 
charge in the underlying quark-pair production and the 
rates of cascade semileptonic decays of 6-hadrons and c- 
hadrons. In regions with no HFR events, we estimate the 
size of this background by extrapolating the HF predic- 
tion from neighboring dimuon control regions that con- 
tain sufficient numbers of events. 

The trilepton HF background is estimated by requiring 
that the normalized HFR sample has a third lepton. If 
there are no events satisfying this requirement, then we 
extrapolate from either neighboring dimuon or trilepton 
control regions with sufficient statistics. For example, 
we have no HFR data in the trilepton signal region. We 
estimate the HF background there by extrapolating from 
the low-^T region, where we have trilepton HFR events. 
For the extrapolation we use the dimuon Ipx distribution, 
using the fact the ]f^T distribution is similar for dimuon 
and trilepton events. We verify this fact with the use of 
MC-simulated bb/cc events. 

For the determination of the systematic uncertainty as- 
sociated with the HF-estimation method, we re-estimate 
the HF background by redefining the HFR dataset using 
either a requirement on the number of silicon detector 
hits for the muon that has large impact parameter (at 
least two silicon hits), and/or applying a requirement on 
the impact parameter significance (|do|/'^Mo| > 5). These 
cuts favor HF events but reduce our HFR dataset statis- 
tics. The HF-estimation method systematic uncertainty 
is about 25% in the signal region. 

Although the HF normalization is extracted from the 
inclusive analysis sample, with dimuon mass greater than 
10.5 GeV/c^ (to avoid the T resonances) and no addi- 
tional Iflx or jet multiplicity requirements, the agreement 
of our HF predictions in both event yields and kinematic 
distributions for all our dimuon and trilepton control re- 
gions is excellent, as we show in Section rVIIl 



V. SUSY SIGNAL SCENARIOS 

This analysis is a generic search for trilepton events 
in which we focus on minimizing the SM background. 
We nevertheless consider two mSUGRA 37] SUSY signal 
scenarios, "SIGl" and "SIG2" , defined by the value of the 
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FIG. 3: Fit of HF+DY+fakes dimuon mass distribution to the observed data for opposite-charge dimuons. The HF normahza- 
tion is the only free parameter of the fit. The blue (light gray) filled histogram is the IfF and the red (dark gray) is the fakes. 
The thick line represents the total background, which is almost exclusively DY, HF, and fakes at the opposite-charge dimuon 
level. The hatched areas indicate the total uncertainty on the prediction. 




FIG. 4: Fit of HF-f DY-|-fakes dimuon mass distribution to the observed data for same-charge dimuons. The HF normalization 
is the only free parameter of the fit. The blue (light gray) filled histogram is the HF and the red (dark gray) is the fakes. The 
thick line represents the total background, which is constituted almost exclusively by HF and fakes for same-charge dimuon 
pairs. The hatched areas indicate the total uncertainty on the prediction. 
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common sfermion mass (mg) and common gaugino mass 
(mi) at unification scale, the trilinear coupling (Aq), the 
ratio of the two Higgs fields vacuum expectation values 
(tan/3), and the sign of the higgsino mixing parameter 
(sign(/i)): 

• SIGl: mo=100 GeV/c^, mi=180 GeV/c^, Ao=0, 
tan/3=5, ^ > 0. The expected cross section a{pp 
xtx2) times the branching ratio B to leptons is 
a X B = 0.642 x 0.22 pb. The cross section was 
obtained using the next-to-lcading order calcula- 
tion of PROSPINO [11] and the branching ratio us- 
ing PYTHIA. The corresponding chargino and light- 
est neutralino masses would be 116 GeV/c^ and 65 
GeV/c^, respectively. 

• SIG2: TOo=74 GeV/c^, mi=168 GeV/c^, Ao=0, 
tan f3— 3, /i > 0. The expected cross section times 
the branching ratio to leptons is axB ^ 1.023 x 0.5 
pb, as given by prospino and pythia. The cor- 
responding chargino and lightest neutralino masses 
would be 103 GeV/c^ and 57 GeV/c^, respectively. 

These two signal scenarios serve as benchmarks of pos- 
sible SUSY signal and were used for the optimization of 
the minimum req uirement in the signal region, which 
is set at 15 GeV [3^. The mass spectrum of the super- 
symmetric particles was obtained with isajet ^3] and 
the events are generated with pythia. The SIGl sce- 
nario leads to three-body decays (Fig. 2a,b) of the light- 
est chargino (xf ) and the next-to-lightest neutralino {X2) 
with branching ratios to electrons and muons suppressed, 
due to the low branching ratio of the gauge bosons to lep- 
tons. On the other hand, the SIG2 scenario leads exclu- 
sively to two-body decays (Fig. 2c, d) of both gauginos to 
sleptons, with xf decaying always to a final-state tau lep- 
ton (produced from a stau decay). Our analysis is more 
sensitive to SIG2, due to the higher cross section and our 
ability to select events with low momentum final-state 
leptons, originating from tau decays. 

VI. SYSTEMATIC UNCERTAINTIES 

The sensitivity of our search to signals of new physics 
and the significance of a potential excess of events are 
influenced by the uncertainties on our background esti- 
mates. Because we perform a counting experiment, we 
concentrate on the uncertainties on the expected number 
of background and signal events. The event-yield system- 
atic uncertainty is naturally different for MC-simulated 
and CDF-data-estimated physical processes. We first 
discuss the systematic uncertainty on the MC-estimated 
backgrounds and SUSY signals and then we treat the sys- 
tematic uncertainty on the CDF-data-based background 
from HF and fake leptons. 

The sources of systematic uncertainty in the signal re- 
gion, with their effect on signal and MC-estimated back- 
ground event yields are: 



• the luminosity uncertainty (6%) [4ll. |43|. 

• the lepton-identification scale factors uncertainty 
(~ 10%), 

• the trigger efficiency uncertainty ('^ 1%), 

• the jet-energy scale uncertainty (~ 1%) [27]; this 
source of systematic uncertainty is responsible for 
migrating events from one control or signal region 
to another, since variations in jet energies affect 
both the corrections to the and the jet multi- 
plicity, 

• the PDF uncertainty (l%-2%) 

• the uncertainty from the theoretical cross-sections 
estimates (5-12% depending on the process) [H, 

• the uncertainty on the initial- and final-state QCD- 
induced radiation (ISR/FSR) j43|, which has an 
effect of 4% and 12% for background and signal 
MC samples, respectively, and 

• the uncertainty induced from the limited MC statis- 
tics: for the SIG2 MC it is ^ 2% for the dimuons 
and ~ 6% for the trileptons; for the standard-model 
background MC it is ~ 3% for the dimuons and 
~ 40% for the trileptons (the latter mainly due to 
the DY-I-7 limited MC statistics). 

All of the above sources of systematic uncertainty are 
correlated among the different physics processes (DY, di- 
boson, tt), with the exceptions of the cross section sys- 
tematic uncertainties and the MC samples' statistical un- 
certainties. Still, the sources of systematic uncertainties 
are uncorrelated with each other and the respective un- 
certainties are summed in quadrature with each other to 
give the total yield uncertainties in the control and signal 
regions. 

The systematic uncertainty associated with the HF- 
estimation method consists of a part that is anti- 
correlated with the DY-|-fakes systematic uncertainty 
(because the HF weights are given by the fit of the 
DY-l-fakes-l-HF to the observed data, and a varied level 
of DY-l-fakes affects these weights) and an uncorrelated 
part (from the fit uncertainty of about 2-4% for the fixed 
DY+fakes level and from the HF-estimation method 
systematic uncertainty, as described in Section FlVB 2p . 
The correlated DY-|-HF-|-fakes systematic uncertainty is 
about 22%. The fake-lepton uncertainty is set to a con- 
servative maximum-envelope 50% level, which is deter- 
mined by studying different jet-triggered CDF samples 
[3^. For the total systematic uncertainty of the pre- 
dicted SM event yield in all control regions we take into 
account all correlations among physics processes. For 
each source of systematic uncertainty affecting the MC 
samples, we vary all MC samples (including DY) in a cor- 
related manner and redo the fit of DY-|-fakes-|-HF to the 
observed data to extract a new HF estimation. The total 
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TABLE I: The dimuon and trilepton event-yield systematic 
uncertainties for the backgrounds in the signal region. The 
uncertainties are summed based on the contributions of the 
separate backgrounds taking into account all correlations. 
The upper part of the table shows the MC-related system- 
atic uncertainties whereas the lower part shows the system- 
atic uncertainties for the CDF-data-estimated backgrounds. 
The uncertainties due to the MC statistics are not shown. 



Source 


Dimuons Trileptons 


Electron scale factors 




±2% 


Muon scale factors 


±8% 


±5% 


Luminosity 


±3% 


±2% 


Trigger efficiency 


±0.5% 


±0.2% 


PDF 


±2% 


±1% 


ISR/FSR 


±2% 


±1% 


Theoretical cross sections 


±3% 


±2% 


Jet-energy scale 


±0.5% 


±0.02% 


Total MC syst. 


±9% 


±6% 


Fakes estimation 


±8% 


±25% 


HF estimation 


±5% 


±2% 


Total (with correlations 


) ±10% 


±24% 



TABLE IL The dimuon and trilepton event-yield systematic 
uncertainties in the signal region, for the SIGl and SIG2 
SUSY scenarios. 



Source 


Dimuons Trileptons 


Electron scale factors 




±6% 


Muon scale factors 


±11% 


±15% 


Luminosity 


±6% 


±6% 


Trigger efficiency 


±0.9% 


±0.5% 


PDF 


±1% 


±1% 


ISR 


±12% 


±12% 


Theoretical cross sections 


±10% 


±10% 


Jet-energy scale 


±0.3% 


±0.6% 


MC statistics 


±2% 


±6% 


Total 


±20% 


±24% 



variation gives us the total effect of the systematic uncer- 
tainty. The same procedure is followed when we include 
the fake-lepton estimation uncertainty and its effect on 
HF due to the above fit. 

The effect of the systematic uncertainties on the back- 
ground and SUSY signal expected event yields in the sig- 
nal region can be found in Tables |T] and |lT] respectively. 



VII. CONTROL REGIONS 

We investigate control regions defined by the dimuon 
mass, IpT, and jet multiplicity, as shown in Fig. [51 Over- 
all, 22 dimuon and trilepton control regions are defined. 
Most control regions we investigate are naturally SM- 
dominated with little expectation of SUSY signal. Low 
IpT and regions are dominated by the HF back- 



ground, whereas the 76 < M^^ < 106 GeV/c^ region is 
almost exclusively populated with Z bosons. The 5 GcV 
gap in the Ifx cuts between the signal and the control 
regions ensures that the low If^T control regions contain a 
negligible amount of signal. We compare the SM event- 
yield predictions with observed events in the control re- 
gions (along with kinematic plots) before looking at the 
signal region. 

We present here a Z boson resonance control region 
("Control z"), a low control region ("Control A"), 
and three control regions ("Control B, C, and D") that 
result from the inversion of one of the three signal region 
requirements at a time [dimuon mass (M^^), or or 
jet multipHcity {N-^cts) respectively]. In region z, we re- 
quire that the muons have opposite charge and that the 
dimuon mass lie between 76 and 106 GeV/c^. We use 
this control region for validating the luminosity, the trig- 
ger efficiencies, and muon-identification scale factors. In 
region A, we require that ]pT < 10 GeV and Af^^ > 10.5 
GeV/c^. This region is used for verifying our knowl- 
edge of HF and fake-lepton backgrounds. In region B, 
we require J^t > 15 GeV, dimuon mass within the Z- 
mass region, and low jet multiplicity (at most one jet). 
This region helps us verify our background prediction in 
a low-yield region as most Z events are characterized by 
low ]^T- In region C, we require < 10 GeV, exclu- 
sion of the Z mass region, M^^ > 15 GeV/c^, and low 
jet multiplicity. This region along with region A are the 
ones with the highest population of HF events. In region 
D, we require > 15 GeV, exclusion of the Z mass re- 
gion, Af^^ > 15 GeV/c^, and high jet multiplicity (more 
than one jet). This region is expected to be the most 
sensitive to tt production. We finally study the dimuon 
events with all signal-region kinematic cuts applied, but 
before the requirement for a third lepton. This is a crit- 
ical control region as the trilepton signal is a subset of 
this region. 

Table Hill shows the expected and observed number of 
dimuon events in our control regions, and Table HVl shows 
the expected and observed number of trilepton events. 
After requiring the presence of a third electron or muon, 
only control regions z, A and C are populated with ex- 
perimental data. Region z trilepton event yields establish 
our understanding of the electron fakes, since the third 
electron in Z boson events is almost exclusively a non- 
prompt electron. On the other hand, trilepton regions 
z and C confirm our understanding of the HF and fake 
backgrounds for trileptons. Figs. El and [71 show the 
dimuon mass and I^t distributions for the dimuon con- 
trol regions. The agreement between observed data and 
prediction in the control regions is satisfactory, both in 
event yields and kinematic distributions. 



VIII. SIGNAL REGION RESULT 

After observing satisfactory agreement between experi- 
mental data and SM predictions in both the dimuon and 
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FIG. 5: The control and signal regions used in our analysis are defined in the dimuon mass vs. I^t plane, with the extra 
requirement of low (< 1) or high (> 1) jet multiplicity. In this paper we show results for the control regions that result from 
the inversion of one of the three main kinematic selections (dimuon mass, missing transverse energy, and jet multiplicity), with 
the addition of a Z-mass control region (Control z) and a 1ow-_^t control region (Control a). The control regions above are 
defined for low jet multiplicity, unless otherwise stated. 



TABLE III: Expected and observed dimuon event yields, in all control regions and the signal region. The expected SUSY signal 
event yield is for the SIG2 mSUGRA scenario. Combined statistical and systematic uncertainties are shown and correlations 
among sources of systematic uncertainty are included. The signal region without a requirement for a third lepton is a dimuon 
control region. 



Region DY HF Fakes Diboson tt Total SM expected SUSY expected Observed 



Control z 


6419 ± 709 




10± 11 


2.4 ±0.2 


1.18 ±0.14 


6433 ± 712 


0.30 ±0.07 


6347 


Control A 


14820 ± 2242 


9344 ± 1612 


2294 ± 1148 


1.03 ± 0.09 


0.12 ±0.03 


26459 ± 1429 


0.9 ±0.2 


26295 


Control B 


217 ±25 




9 ± 7 


1.7 ±0.2 


0.27 ±0.05 


227 ± 26 


0.5 ±0.1 


253 


Control C 


5770 ± 1043 


2238 ± 384 


466 ± 234 


0.49 ± 0.07 


0.02 ±0.01 


8474 ± 857 


0.7 ±0.2 


8205 


Control D 


7.8 ± 1.5 


9±4 


0.3 ±0.3 


0.21 ±0.07 


4.1 ±0.4 


22 ± 5 


1.8 ±0.4 


23 


Signal Reg. 


169 ± 30 


90 ±20 


49 ± 25 


6.5 ±0.4 


0.96 ±0.11 


315 ± 37 


17 ±3 


297 



TABLE IV: Expected and observed trilepton event yields, in all control regions and the signal region. The expected SUSY signal 
event yield is for the SIG2 mSUGRA scenario. Combined statistical and systematic uncertainties are shown and correlations 
among sources of systematic uncertainty are included. 



Region 


DY 


HF 


Fakes 


Diboson 


tt 


Total SM expected SUSY expected Observed 


Control z 


0.2 ±0.2 




2.5 ± 1.2 


0.26 ± 0.06 




3±1 


0.06 ±0.01 


4 


Control A 


0.3 ±0.2 


6 ± 3 


7.6 ±3.8 


0.25 ± 0.08 




14 ±4 


0.08 ± 0.02 


16 


Control B 






0.2 ±0.1 


0.094 ± 0.009 




0.3 ±0.1 


0.10 ±0.03 





Control C 


0.2 ±0.2 


3 ± 2 


2±1 


0.10 ±0.06 




5 ± 2 


0.06 ± 0.02 


8 


Control D 






0.02 ± 0.01 


0.003 ± 0.002 


0.011 ±0.008 


0.03 ± 0.01 


0.04 ±0.02 





Signal Reg. 




0.06 ±0.04 


0.2 ±0.1 


0.15 ±0.06 




0.4 ±0.1 


1.7 ± 0.4 


1 



trilepton control regions, we look at the CDF data in 
the trilepton signal region. We observe one event con- 
taining three muons (trimuon event). The event is char- 
acterized by low track activity and three well-identified 
muons that are produced within 40 degrees in 4>. Two 
of the muons are energetic, with transverse momenta of 
45 and 21 GeV/c, and the third one is a soft muon with 



Pt of 8 GeV/c. Table fVl shows the main properties of 
this event. It is interesting to note the close values of all 
three dimuon masses. The event includes one hadronic 
jet and two energy clusters of mostly electromagnetic en- 
ergy with transverse energy of ~ 41, ~ 9, and ~ 4 GeV 
respectively. The jet and the muons originate from the 
same and only high-quality primary vertex. If the elec- 
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FIG. 6: Dimuon mass and distributions for the SM background in the dimuon control regions z (a,b), A (c,d), and B (e,f). 
The background histograms are stacked. The CDF data are indicated by points with error bars. 



tromagnetic energy clusters correspond to real photons, 
the event would also be interesting in the gauge-mediated 
supersymmetry breaking (GMSB [44]) scenario, where 
the lightest neutralino decays to a photon and a grav- 
itino, which is the LSP. In that case, the final leptonic 
signature of the chargino- neutralino production would be 
three leptons, two photons, and ]^t- 

Fig. [5] shows where the one trimuon event observed in 
the signal region appears in the expected distributions of 
kinematic variables for the signal and the backgrounds. 
Kinematic distributions include the three-dimensional 
opening angle between the leading muons, Aip. Figs. 
1^ and [TUl show the transverse and lego detector displays. 



respectively, for this trimuon event. The Poisson proba- 
bility to see one event or more, when we expect 0.4±0.1, 
is 32.6%. 

It is interesting also to interpret this event in the con- 
text of a search only for trimuon events. The diboson 
backgrounds remain, but the large source of fakes in the 
dimuon-|-e sample is reduced. The total trimuon back- 
ground estimation is 0.16 ± 0.04 events. The Poisson 
probability to observe one event or more, when we ex- 
pect 0.16 ± 0.04, is 14.7%. We conclude that our event 
yield is statistically consistent with the SM prediction, 
noting that most of the kinematics of the event - espe- 
cially the three-dimensional opening angle of the leading 
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FIG. 7: Dimuon mass and for the SM background in the dimuon control regions C (a,b), D (c,d), and dimuon signal region 
(e,f). The background histograms are stacked. The CDF data are indicated by points with error bars. 

TABLE V: Observed trimuon event properties. 

Kind of muons CMUP - CMX - CMX 

Pt of muons (GeV/c) 45.0, 21.1, 7.8 

■q of muons -0.2, -0.9, 0.8 

<j> of muons (deg.) 359, 321, 340 

Isolation of muons (GeV) 2.4, 0.2, 1.1 

Charge of muons -1, 1, -1 
Dimuon masses (GeV/c^) 29.3(1&2), 21.7(1&3), 25.7{2&3) 

Transverse mass (muon-f.^y) 86.4, 51.4, 34.2 

3-d Ai(5(leading muons) (deg.) 46.3 

$T{GeV) 43.8 

$T 4> (deg.) 205.6 

Number of Jets 1 

of jet (GeV) 41.1 

rj of jet -1.6 

<j> of jet (deg.) 102.9 
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FIG. 8; Kinematic variables for the SM background and the SIG2 SUSY signal, in the trilepton signal region. The background 
histograms are stacked; the signal histogram is not. The CDF data are indicated by points with error bars. 
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FIG. 9: The trimuon event in the transverse view of the cen- 
tral CDF detector. Tracks with transverse momenta above 
1 GeV/c are shown. 

muons Ai^ - are consistent with new physics expectation, 
as can be seen in Fig. [S] 

We have combined the results of this analysis with 
other CDF trilepton analyses to set exclusion limits in 
several models. For mSUGRA with no slepton mixing, 
we set a lower limit for the chargino mass of 129 GeV/c^, 
which corresponds to an upper limit in a x S of about 
0.25 pb at the 95% confidence level [13]. 
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